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We analyze the DC charge transport in the quantum critical regime near a d-wave 
Pomeranchuk instability in two dimensions. The transport decay rate is linear in 
temperature everywhere on the Fermi surface except at cold spots on the Brillouin 
zone diagonal. For pure systems, this leads to a DC resistivity proportional to T'^/^ 
in the low-temperature limit. In the presence of impurities the residual impurity 
resistance at T = is approached linearly at low temperatures. 
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Pomeranchuk instabilities Jj leading to symmetry-breaking deformations of the Fermi 
surface in interacting electron systems have attracted much interest in the last few years. 
Interactions favoring a Pomeranchuk instability with d-wave symmetry have been found in 
the two most intensively studied single-band models for cuprate superconductors, that is. 



the two-dimensional 



t-jQ 



and Hubbard 



^ model. These models thus exhibit enhanced 
"nematic" correlations, as usually discussed in the context of fluctuating stripe order 
5j. Signatures for incipient nematic order with d-wave symmetry have been observed 
in various cuprate materials 6|]. In particular, nematic correlations close to a d-wave 
Pomeranchuk instability provide a natural explanation for the relatively strong in-plane 
anisotropy observed in the magnetic excitation spectrum of YBa2Cu30y j3, IJ. A spin 
dependent Pomeranchuk instability was recently invoked to explain a new phase observed 
in ultrapure crystals of the layered ruthenate metal Sr3Ru207 9|, and also to account for 



12 Q. 



ity provide an interesting route to 
The properties of single-particle 



a puzzling phase transition in URu2Si 

Critical fluctuations near a Pomeranchuk instabi 
non- Fermi liquid behavior in two dimensions jlll . [l2 |. 
excitations near a quantum critical point associated with a Pomeranchuk instability have 
been studied already in considerable detail jisl. [l^. For a d-wave Pomeranchuk instability 
in an electron system on a square lattice the singular part of the electronic self-energy 
is proportional to (i^, where d]^ is a form factor with d-wave symmetry [12]. At the 
quantum critical point, the real and imaginary parts of the self-energy scale as \uj\'^^^ 



with energy 



11 



12|. This leads to a complete destruction of quasi-particles near the 
Fermi surface except for the "cold spots" on the Brillouin zone diagonal, where the form 
factor dk vanishes. In the quantum critical regime at T > the self-energy consists of 
a "classical" and a "quantum" part with very different dependences on T and uj. The 
classical part, which is due to classical fluctuations, dominates at = and yields a 
contribution proportional to \/T/ logT to the imaginary part of the self-energy on the 



Fermi surface 13|. 



In this letter we compute the temperature dependence of the DC resistivity in the 
quantum critical regime near a Pomeranchuk instability in two dimensions. We obtain a 
momentum dependent transport decay rate J^iT) which is linear in temperature for all 
momenta on the Fermi surface except at the cold spots on the Brillouin zone diagonal. 
Adding a conventional T^-term to 7^ (T) we obtain an overall resistivity p{T) proportional 
to T'^/^ at low temperatures. In the presence of impurities, the residual resistivity at zero 
temperature is approached linearly. 



Our calculations are based on a phenomenological lattice model |l2[ |. 



H = Hq + ^ /kk'(q) nk(q) nk'(-q) , (1) 



k,k',q 

J 



where Hq is a kinetic energy, n]^{q) = '^k-q/2 a'^'k+q/2,a, and V is the volume of the sys- 
tem. Since the Pomeranchuk instability is driven by interactions with small momentum 
transfers (forward scattering), we choose a coupling function /kk'(q) which contributes 
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only for relatively small momenta q. This excludes other instabilities such as supercon- 
ductivity or density waves. We consider an interaction of the form Q 

/kk'(q) = + ^(q) 4 4' (2) 

with u{q) > and g{q) < 0, and a form factor with d^'^-y'^ symmetry, such as (ik = 
cosA;^. — cos ky. The coupling functions u{q) and g{q) vanish if |q| exceeds a certain 
small momentum cutoff Or- This ansatz mimics the effective interaction in the forward 
scattering channel as obtained from renormalization group calculations [3[ for the two- 
dimensional Hubbard model. The uniform term originates directly from the repulsion 
between electrons and suppresses the electronic compressibility of the system. The d- 
wave term drives the Pomeranchuk instability. 

The Pomeranchuk instability can actually be preempted by a first order transition 
at low temperatures, where the Fermi surface symmetry changes abruptly before the 



fluctuations become truely critical 



161 ]. However, for reasonable choices of hopping and 



interaction parameters the system is nevertheless characterized by strong fluctuations on 



the symmetric side of the transition 15|. The first order character of the transition is 
suppressed by the uniform repulsion u in (2), and for a favorable but not unphysical choice 



15| 



of model parameters a genuine quantum critical point can be realized 

Near the Pomeranchuk instability, the electrons interact via a singular effective inter- 

nn 

action of the form [12, 113II 



gdkd], 

(WO' + ^o'|qP-^K(c|q| 



rkk'(q, = ,^ , ^sui-^ ,r..//„i„iM ) (3) 



where q and v is the momentum and energy transfer, respectively. The parameters 
g = g{0), ^0 and c can be treated as constants, whereas the correlation length ^ depends 
sensitively on control parameters and temperature. In the quantum critical regime C,{T) 
is proportional to (T| logT|)~^/^. 



The electron self-energy S(k, u;) has been computed previously 12, llSl, in random 



phase approximation (RPA) with the effective interaction (3). In the quantum critical 
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regime one finds 



13| 



ImE(kp,0) = ^^Te(T) (4) 

for momenta on tlie Fermi surface (fk^ is the Fermi velocity). Tlie corresponding 
approximation for the electrical resistivity involves the RPA self-energy and current ver- 
tex corrections due to particle-hole ladder diagrams, in close analogy to the Born ap- 
proximation for impurity scattering [l7|. We assume that the Pomeranchuk fluctuations 
thermalize sufficiently rapidly such that the effective interaction (3) is not modified by 
the electric current. This relaxation to equilibrium is not described by the model (1) and 
las to be provided by additional terms such as umklapp, impurity or phonon scattering 



18|, 



19|. 



The DC conductivity (inverse resistivity) can be obtained from the retarded current- 
current correlation function 11 as 

a^i = — lim lim — lmlljii(q,uj) (5) 

For cubic symmetry the conductivity tensor is diagonal. The current-current correlator 
can be expressed (exactly) in terms of single-particle Green functions G and current 
vertices. Performing an analytic continuation from Matsubara to real frequencies and 
taking the DC limit, one obtains 



~ J dujf'iu;) J ^ A;(k) |G(k,c^)|2A,,(k,c.) , (6) 



where f{uj) is the Fermi function, A°(k) = = Vek the bare current vertex, and 
A(k, cu) is the interacting current vertex in the mixed advanced-retarded DC limit, that 
is, A(k,co') = A(k,ci; + iO^; ]<., uj — iO^) . The product IGp can be expressed in terms of the 
single-particle spectral function A(k,u) = — ;ilmG'(k, c<j) and the (retarded) self-energy 

as 



At low temperatures the derivative of the Fermi function f'{uj) has a sharp peak of width 
T at = 0. Since all other factors under the integral in Eq. (6) have a broader uj- 
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dependence, one can replace f'{uj) by —S{uj), such that the conductivity can be written 

as 

'---^i^.^>)^^^AKO). (8) 

The interacting current vertex includes all particle-hole ladder vertex corrections. It is 
thus obtained from a linear integral equation, which can be written as 

72. 



X Imrkk(q, 6) -4(k + q^^ + ^) A(k + q, c + e) (9) 
lmL(k + q, + e) 

after analytic continuation to the real frequency axis. Here 6(e) is the Bose function. 

For T the correlation length ^(T) diverges. Repeating the arguments used for 
the calculation of S(k, tu) in Ref. [l3|, one finds that the integration variable e in Eq. (9) 
scales as and can therefore be set to zero in the arguments of A, E, and A on the 
right hand side of Eq. (9). Expanding the Bose function as b{e) ~ T/e, one can carry out 
the e-integration explicitly, / c/e Imrkk(q, e) = vrrkk(q, 0), yielding a closed equation 
for the static current vertex A(k) = A(k, 0) 

72, 



(10) 



The same result could have been obtained by considering only the classical fluctuations, 
that is, by including only the term rkk(q, "^en) with Matsubara frequency e„, = in the 
Matsubara sums for the current vertex corrections. At this point the equations for ajj/ and 
A(k) are formally identical to those obtained from the Born approximation in disordered 
electron systems with a k-dependent long-ranged disorder correlator given by rkk(q, 0). 

Inserting the ansatz A(k) = A(k)vk into the equation for the current vertex, one 
obtains the following equation for the function A(k): 

A(k) = 1 + r / ^r,,(q. 0) °> ■ A(k + q) . (11) 

J [271)^ ImS(k + q, 0) 

Since the conductivity is dominated by momenta near the Fermi surface, we now focus 
on the case k = k^. For large ^ the above integral is dominated by small momentum 
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transfers q of order due to the effective interaction Fkklq, 0). The spectral function is 
peaked for momenta on the Fermi surface, with a width determined by ImS(k^, 0), which 
is proportional to T^(T). The self-energy S(k, 0) varies on a momentum scale of order 
for momentum shifts perpendicular to the Fermi surface |13i]. The same can be expected 
for A(k), since the current vertex correction can be related to the shift of the self-energy 
in the presence of a field coupled to the current operator. Since T^^(T) oc 1/logT in 
the quantum critical regime, and since the tangential q-dependence of ImS(kp' + q) and 
A(ki? + q) is negligible on the scale we may neglect the q-dependence of lmI](kF + q) 
and A(ki? -|- q) in (11) altogether, which can then be solved explicitly, yielding 



A(k^) 
Using 



1 - 



vrT 



ImS(ki7,0) J (27r) 



'^''^ Fk,k, (q, 0) A{kF + q, 0) ""^^ ' ^^-^^ 



/d 
^Fk^k^(q, 0)A(kf + q,0) , 

which is true within self-consistent RPA restricted to classical fluctuations 
write Akf, as 

A(k^) = 7k,/7k; , 
where ■jkp = — ImS(kir, 0) is the single-particle decay rate while 



(12) 



(13) 



one can 



-vrT 



Fk^k^(q,0)A(k^ + q,0) 1- 



Vkf. ■ Vk^+q 



(14) 



(15) 



is the scattering rate relevant for transport. 

The momentum integral in the expression (8) for the conductivity is peaked at the 
Fermi surface. For T with ^(T) oc (TlogT)^^/^ one can replace 74(k, 0) under 
the integral by (5(ek — /i), neglecting possible corrections of order 1/logT, such that the 
conductivity can be written as a Fermi surface integral. Inserting Eq. (14) for Xlkp), one 
obtains 



87r2 



■F tr 

7k. 



(16) 



for the diagonal part a = ajj of the conductivity tensor. 
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To compute 7^''^, we parametrize the (small) momentum transfer q in Eq. (15) by 
radial and tangential components, g,. and g^, respectively. For T 0, we may again 
approximate ACkp + q, 0) by a (5-function, 5(ek^_|_q — fj,). The dispersion relation can be 
expanded as Ck^+q - /i = i^kF^'r + Qt / ('^''^kp) ^ where (m^^)"^ = (9|^ek|,^^ . Since kp + q 
is confined to the Fermi surface, the momentum transfers q are predominantly tangential 
to the Fermi surface in hp, such that the term of order cannot be neglected compared 
to the term linear in g^. After expanding also the kinematic factor 1 — ^'^pJ'^F+'i linear 
order in and quadratic order in gj, the momentum integral in Eq. (15) can be performed 
analytically. For ^~^(T) ^ ^kp'Vkp, one obtains 

7t', = J|<arctau(^-^)/<-,,4,r. (17) 

where qc is the momentum cutoff and 

1 /^Vk ■ Sfc^Vk „2 



= TTi: 7 Vk • Sfc^Vk 

2vi \ Vkml ' y ^^^^ 

The function K^^^ has units of inverse momentum squared. For a quadratic dispersion 

relation, = /c^/(2m), one has = m and K\f_p = l/(2/c|,). The scattering rate 7k^ 

is thus linear in T at low temperatures. Note that the correlation length ^ (T) does not 

appear in the asymptotic low temperature behavior of 7k^. This behavior in the quantum 

critical regime can be contrasted with the behavior in the Fermi liquid regime close to 

the quantum critical point. For the latter case 7k^ is proportional to ^^T^logT, where 

e=e(T^o) y. 

Due to the d-wave form factor the prefactor of the T-linear behavior of 7k^ varies 
strongly along the Fermi surface and vanishes on the Brillouin zone diagonal. This is 
reminiscent of the cold spot scenario of transport in cuprates 20(]. Inserting 7^"" from (17) 



in Eq. (16) for the conductivity one obtains a divergent Fermi surface integral, due to 
the zeros of 7k^ at the cold spots k^. In the absence of any other scattering mechanism, 
the conductivity would thus be infinite. However, other (than d-wave forward scattering) 
residual interactions will lead at least to the conventional Fermi liquid decay rate of order 
all over the Fermi surface, including the cold spots. Including a Fermi liquid term 
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of order T^, the scattering rate has the form 'j^^iT) = a^^pT'^ + bi^^dl^^T, where the 
coefficients ai^^ and b^^ are finite for all hp. Inserting this ansatz into Eq. (16), one finds 
a resistivity 

p(T) = ^v{^T3/^ (19) 

for low T. Taking the well-known logarithmic correction to the T^-behavior of the scatter- 
ing rate in two-dimensional Fermi liquids into account, one obtains p{T) oc T^/^ | logT|^/^ . 

In the presence of impurities, the scattering rate has the form 7k^(T) = Tk^'^ + feicF'^k^T", 
for temperatures low enough that the Fermi liquid term of order can be neglected 
compared to the impurity term. For T ^ one then obtains a finite residual resistivity 
determined exclusively by impurity scattering. For low finite temperatures the resistivity 
increases linearly with T as long as T ^ Ik^^/bkp. For T ^ Y^/bkp one obtains 
p(T) oc T^/^, with a prefactor proportional to (7^™^ &kj,)"'^^^/'i^kj, , provided that impurity 
scattering still dominates over the conventional Fermi liquid contribution to '-/^^{T). 

In summary, we have analyzed the DC charge transport in the quantum critical regime 
near a d-wave Pomeranchuk instability in two dimensions. It turned out that the relax- 
ation of the electric current is dominated by classical fiuctuations. The transport decay 
rate jI^^ (T) is linear in temperature everywhere on the Fermi surface except at cold spots 
on the Brillouin zone diagonal. For pure systems, this leads to a DC resistivity propor- 
tional to T^/^ in the low-temperature limit. It is tempting to associate this result with 
the unusual T^/^-law observed for the resistivity in overdoped La2-xSrxCu04 2l(]. In the 
presence of impurities the residual impurity resistance at T = is approached linearly at 
low temperatures. 
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